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Abstract 

We present an update of the study of the Y(4260) resonance, produced in the process e + e~ — > 
HSR TT + TT^J/ifj using initial-state radiation events at the PEP-II e + e~ storage rings. This study 
is based on 454 fb^ 1 of data recorded with the BABAR detector at a center-of-mass energy in the 
T(4S) resonance region. From a fit with a single non-relativistic Breit-Wigner shape we obtain 
updated parameters for the Y(4260) resonance which are my = 4252 ± 6 ±1 MeV/c 2 and T Y = 
105 ± 18+g MeV/c 2 ; we also measure B(tt + tt- J/tp )T e+e - = (7.5 ± 0.9 ± 0.8) eV. We cannot confirm 
the recent BELLE observation of a broad structure around 4.05 GeV/ c 2 in this decay mode. 



Submitted to the 34 International Conference on High-Energy Physics, ICHEP 08, 
29 July — 5 August 2008, Philadelphia, Pennsylvania. 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



Work supported in part by Department of Energy contract DE-AC02-76SF00515. 



1 INTRODUCTION 



The observation of X(3872) pQ, followed by the discovery of many other states such as the Z(3930) 
[2], the Y(3940) [3j and the X(3940) [3], has reopened interest in charmonium spectroscopy. Some 
of these resonances cannot be fully explained by a simple charmonia model [5]; four-quark state [6] 
and D°D°* molecule [7] are some of the interpretations that have been proposed to explain their 
nature. Among these new states the Y(4260) [8j that has been observed by BABAR in the process 
e + e~ — > 'Jisr Y(4260)— > ^isr J/vpir + ir~ , where ISR denotes initial state radiation. Being formed 
directly in e + e~ annihilation, it is known to have J PC = 1 . Nonetheless its properties do not fit 
with any simple charmonium interpretation and its nature is still unclear. A recent ISR analysis 
by BELLE [9] suggests the existence of a second state in the same decay mode, called Y(4008) [10J, 
with the same production mechanism. 

After the discovery of Y(4260), other J PC = 1 states, produced through the same e + e~ 
initial state radiation mechanism, have been observed in the ip(2S)ir + ir~ final state, such as the 
Y(4320) and the Y(4660) PHH2J. 

In this letter, we present udated measurements of the resonance parameters and B(it + it~ J/ip )T e + e - 
of the Y(4260) based on the full BABAR data set which is approximately two times larger than the 
one used in the original analysis [8]. 

2 THE BABAR DETECTOR AND DATASET 

The original BABAR analysis [8] is updated here with the same method on a data set corresponding 
to a total integrated luminosity of 454 fb -1 obtained at the SLAC PEP-II B-factory, running at the 
center-of-mass energy near the T(4S) resonance. 

The BABAR detector is described in detail elsewhere |13| . Charged-particle momenta are mea- 
sured in a tracking system consisting of a five-layer double-sided silicon vertex tracker (SVT) and a 
40-layer central drift chamber (DCH), both situated in a 1.5-T axial magnetic field. An internally 
reflecting ring-imaging Cherenkov detector (DIRC) provides charged-particle identification together 
with dE/dx measurements from SVT and DCH. A CsI(Tl) electromagnetic calorimeter (EMC) is 
used to detect and identify photons and electrons. Muons are identified using information from the 
Instrumented Flux Return (IFR) system, together with E/p, where the energy E is determined by 
the EMC and the momentum p by the SVT and DCH. 

3 ANALYSIS METHOD 

In this analysis we reconstruct events: 

e + e~ -> 7 ISR 7r + 7r~ J/tp, (1) 

where 7 ISE represents a photon that is radiated from the initial state , lowering the center-of-mass 
energy of the e + e~ system. This mechanism allows to perform an energy scan of e + e~ interactions 
at the B-factories, where the center-of-mass energy of the interactions is fixed in the region of the 
T(4S) resonance. 

A candidate J/tjj meson is reconstructed via its decay to e + e~ or The lepton tracks, at 

least one of which must be identified as an electron or muon candidate, must be well reconstructed 
and must originate from the same vertex: a geometric fit of the J/ip candidate is conducted with 
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beam-spot constraint. An algorithm to recover energy lost to bremsstrahlung is applied to electron 
candidates. An e + e~ pair with its invariant mass within the interval (—75, +40) MeV/c 2 around the 
PDG [13] J/0 mass is taken as a J/0 candidate. For a pair, the interval is (—40, +40) MeV/c 2 . 

The J/0 candidate is then kinematically constrained to the nominal J/0 mass and combined with 
a pair of oppositely-charged tracks identified as pion candidates. 

In the event reconstruction, we do not require observation of the initial state radiation photon 
7 ISR , as it is preferentially produced along the beam directions and it is not detected. 

We select e + e~ — > 7 ISR 7r + 7r~ J/0 events with the following criteria: the invariant mass squared re- 
coiling against the tt + it~ J/0 system (m 2 ec ) is required to be in the range (—0.50, +0.75) ( GeV/c 2 ) 2 ; 
the total number of reconstructed tracks in the event < 5. Since the transverse component of the 
missing momentum {p^pmiss) ls sman f° r events with initial state radiation, we select events with 
PTmiss < 2.25GeV/c. The cosine of the angle between the £ + momentum in the J/0 rest frame 
and the J/0 momentum in the center-of-mass frame (cos 6i) is required to be | cos#/| < 0.925. 

After this selection, a clean ip(2S) signal is apparent in the ir + ir~ J/ip invariant mass spectrum, 
as shown in figure [2 in the same distribution, a clean Y(4260) peak is also visible. 




5 5.2 5.4 

my/xifJrtrXGeVAr) 



Figure 1: The invariant mass distribution of tt + tt J/0 candidates with a logarithmic vertical scale. 

By means of Monte Carlo simulations we determine that the selection efficiency at the Y(4260) 
energy is 15.3%. The method used for the efficiency determination is checked at the '0(25'), where 
we obtain a value of the cross section that is consistent with the measurement performed by previous 
experiments [9j. The ip(2S) mass is measured as 3685.35+0.02 MeV/ c 2 (where the error is statistical 
only); the mass shift with respect to the PDG value [H] is taken into account as a systematic error 
in the Y(4260) mass measurement. 

The Monte Carlo 7r + 7r~ J/0 invariant mass resolution and mass scale have been calibrated by 
comparing the widths of Tr + ir~Jfy invariant mass distributions from ip(2S) decays in data and 
Monte Carlo. We find that the Monte Carlo simulation reproduces the observed r.m.s. of the 
ir + ir~J/ip distribution for the tp(2S) state. The mass resolution is around 5MeV/c 2 in the mass 
range 4.16GeV/c 2 < m(7r+vr- J/0 ) < 4.36GeV/c 2 

The typical ISR-production signature of the Y(4260) can be obtained by subtracting distri- 
butions for events with -k + tt~ J/ip mass in the region (3.95,4.1) GeV/c 2 and (4.4,4.55) GeV/c 2 from 
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Figure 2: The distribution of m^ ec . The points represent the data events passing all selection 
criteria and having a ir + ir~J/ip mass near 4260MeV/c 2 . They are obtained by subtracting the 
distribution from neighboring tt + tt~ J/ip mass regions (see text). The solid histogram represents 
ISR Y(4260) Monte Carlo events. 



those with mass in the signal region defined as (4.1,4.4) GeV/c 2 . The distributions of m 2 ec and of 
the cosine of the 7r + 7r~ J/ip system polar angle in the cm. frame are shown in figure [2] and [3] along 
with the corresponding distributions for ISR Y(4260) Monte Carlo events. 
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Figure 3: The distribution of cos8* +n _ ju . The points represent the data events passing all selection 
criteria and having a ir + ir~J/tp mass near 4260MeV/c 2 . They are obtained by subtracting the 
distribution from neighboring tt + tt~ J/ip mass regions (see text). The solid histogram represents 
ISR Y(4260) Monte Carlo events. 
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4 RESULTS 



We have performed an unbinned maximum likelihood fit to the tt + it~ J/i/j invariant mass distribution 
between 3.8GeV/c 2 and 5GeV/c 2 . The signal probability density function (PDF) is described by 
convolving a non-relativistic Breit-Wigner function with a Gaussian and we use a first degree 
polynomial function to describe the background. All the fit parameters for the signal and the 
background PDF are floating, except the Gaussian sigma which is fixed to the resolution found at 
the i^{2S) in data but linearly scaled to the Y mass value based on the Monte Carlo simulation. 

From the fit shown in figure [4] we obtain a signal yield Ny = 344 ± 39 and resonance mass 
m Y = 4252 ± 6 t\ MeV/c 2 and width T Y = 105 ± 18 tgMeV/c 2 . Systematic uncertainties include 
contributions from the fitting procedure (evaluated by changing the fit range and the background 
PDF), the mass scale, the mass-resolution function and the dependence on the model of Y(4260)— > 
ir + iT~J/ip decay. They have been added in quadrature. 

To measure ^(vr+vr- J/ip )T e+e - we use Eq. El where N(j^(2S)), NfrY), m(ip(2S)), m(Y), 
e(ijj(2S)), e(Y) and W(ip(2S)), W(Y) are the number of events, mass, selection efficiency and the 
photon emission probability density function for the ip(2S) and Y(4260) respectively. 

r ee (y)g(y^7r+7r-JA/Q 

T ee ($(2S))B(${2S)^ TT+ir-J/i;) 

( NjjY) \ / m(Y) x ( eW2S)) \ , W(i/,(2S)) \ 
VJV( 7 ^(25))/ \m(il](2S))) V s(Y) J V W(Y) J { ' 

This way, the entire uncertainties of integrated luminosity and B(J/?p — > are canceled. Most 

uncertainties of the selection efficiency, particle ID efficiency, tracking efficiency, and photon emis- 
sion probability density W(s,x) are canceled out. Meanwhile, we also introduce some new uncer- 
tainties pertaining to the ISR tp(2S) such as B(ip(2S) — ► ir + ir~ J/ip), T^^sy^e-e- > an d statistical 
uncertainty of the N(ip(2S)). We obtain B(tt + tt- J/i/j)T e+e - = (7.5 ± 0.9 ± 0.8) eV. 




m(7i + 7iJ/^)(GeV/c 2 ) 

Figure 4: The -k + tt~J/iJ) invariant mass spectrum in the range 3.8-5.5 GeV/c 2 . The points with 
error bars represent the selected data, the solid curve shows the result of the fit described in the 
text, while the dashed curve represents the background component. 
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A recent BELLE analysis reported evidence for a broad enhancement with mass m = 4008 ± 
MeV/c 2 and width T = 226 ± 44 ± 87MeV/c 2 [5]. We cannot confirm this observation and 
obtain an upper limit B(ir + ir^ J/tp)T e + e - < 0.7 eV at 90% C.L. for this state. 

We also study the tt + tt~ invariant mass distribution for Y(4260) events, obtained by dividing 
the total sample into several dipion invariant mass regions and fitting the Y(4260) peak in each 
one of these regions. 

The results reported in figure [5] show that, in Y(4260) decays, the tt + tt~ invariant mass distri- 
bution does not agree with the MC based on phase space and tends to peak at large values. 
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Figure 5: Di-pion invariant mass distribution of Y(4260) events: the points with error bars are 
signal events, the histogram is the phase-space distribution from MC events. 



5 CONCLUSIONS 

In summary, we have analyzed initial-state radiation events to study the process e + e~ — * Tr + ir~J/'tjj 
across the charmonium mass range. We observe 344 ± 39 Y(4260) events with my = 4252 ± 

6 i 2 MeV/c 2 and T Y = 105 ± 18 tgMeV/c 2 . The ir + ir~ invariant mass distribution, in Y(4260) 
region, tends to peak at large values consistently with other studies [8] [9]. There is no evidence for 
the broad enhancement reported by BELLE around 4.05GeV/c 2 [9] and we obtain an upper limit 
S(vr + 7r-J/^)r e+e - < 0.7 eV at 90% C.L. for this state. 
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